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in  the  limit  as  the  small  positive  parameter  t -*■  0 . Any  acquaintance  with 
singular  perturbations  would  lead  one  to  expect  that  the  unique  solution 
of  this  initial  value  problem  might  converge  as  € -*■  0 to  a solution  Zq 
of  the  limiting  equation 

(3)  f(ZQ,t,0)  = 0 

away  from  any  "boundary  layer"  regions  of  nonuniform  convergence..  In  the 

singular  situation  that  f (Z  ,t,0)  is  a singular  matrix,  this  limiting 

z u 

equation  has  an  infinite  number  of  solutions,  sc  the  limiting  equation  (3) 
is  not  adequate  for  the  purpose  of  determining  the  reduced  problem  satis- 
fied by  the  limiting  solution  (assuming  that  such  a limit  exists).  Tradi- 
tional methods  apply  whenever  f is  everywhere  stable.  The  solution  to 

z 
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(l)-(2)  then  converges  to  the  unique  solution  of  (3)  av/ay  from  the  initial 
poirit  t = 0 (cf . , e.g.,  O'Malley  (197*0).  Indeed,  the  solution  can  then 
be  uniformly  asymptotically  represented  as  the  sum  of  a function  of  t and 
a decaying  function  of  the  stretched  variable  t = t/£  . 

Many  singular  perturbation  problems  can  be  written  in  the  form  (l)-(2). 
Indeed,  partial  differential  equations  are  often  disguised  as  such  ordinary 
differential  equations  in  appropriate  Banach  spaces  (cf.  Friedman  (1969)) 
and  boundary  value  problems  for  ordinary  differential  equations  are  often 
conveniently  studied  via  shooting  techniques  (i.e.,  by  using  solutions  of 
auxiliary  initial  value  problems)  (cf.  Bernfeld  and  Lakshr.ikantham  (197*0). 
Singular  problems  (i.e.,  those  where  f is  singular)  occur  throughout 
the  literature,  though  they  do  not  seem  to  have  been  analyzed  in  a system- 
atic fashion.  We  note  that  Fife  (1972)  studied  problems  "whose  degenerate 
forms  have  many  solutions".  His  work  was  motivated  by  the  need  to  analyze 
the  dispersion  of  flow  in  small  tubes  (cf.  Fife  and  Nicholes  (1975)). 
Likewise,  Gordon  (1975)  considered  the  quasilinear  problem  (1)  in  a Banach 

space  where  the  matrix  f (z,t,0)  = f (0,t,0)  had  a one-dimensional  null 

z z 

space.  His  applications  include  a nonlinear  stability  analysis  of  the  con- 
vection state  in  a Benard  problem  (cf.,  Gordon  and  Hoppensteadt  (19  )). 

Our  interest  was  generated  by  studying  singular  arc  problems  in  optimal 
control  (cf.  O’Malley  (1976)).  The  solution  of  these  singular  singular- 
perturbation  problems  involves  projections  and  a resulting  reduction  to 
certain  algebraic  systems  and  lower  dimensional  systems  of  differential 
equations,  both  for  the  outer  solution  (asymptotically  valid  away  from 
regions  of  nonuniform  convergence)  and  for  the  boundary  layer  corrections 
(which  allow  nonuniform  convergence).  Such  considerations  are  basic  to  our 
general  plan  of  attack  and  to  the  related  approach  of  Butuzov  and  Vasil 'eva 
(1970)  and  Vasil'eva  (1975). 

We  shall  consider  the  initial  value  problem  (l)-(2)  under  the  primary 
assumptions  that  the  matrix 

(*0  fz(z»t,0)  = fQz(z,t) 

has  constant  rank  k,  0 < k < M,  that  its  null  space  is  spanned  by  M - k 

linearly  independent  eigenvectors,  and  that  f has  k stable  eigenvalues 

z 

(counting  multiplicities)  for  all  z and  t.  Assuming  sufficient  dif- 
ferentiability of  f and  asymptotic  expansions  for  both  f and  z^,  we 
shall  attempt  to  construct  an  asymptotic  solution  of  the  form 


z(t ,c ) = z(t  ,£ ) + n (t  ) 


(5) 


whore  Z and  n have  asymptotic  power  series  expansions  in 
terms  of  II  all  tend  to  zero  as  the  stretched  variable 


and  the 


(6) 


t = t/£ 


tends  to  infinity.  Since  Z(t,£)  provides  the  asymptotic  solution  for 
t > 0,  it  is  called  the  outer  solution,  while  II  (x, 6)  is  called  the 
boundary  layer  correction  since  it  is  asymptotically  negligible  away  from 
the  initial  region  of  nonuniform  convergence.  Simple  examples  show  that 
the  form  of  solution  (5)  would  differ  if  there  were  a deficiency  of  eigen- 
vectors in  the  null  space  of  f^^  (cf.  Butuzov  and  Vasil 'eva  (1970)). 

To  proceed,  it  is  helpful  to  realize  that  the  matrix  f0z(z,t)  can 

be  put  into  its  row-reduced  echelon  form  by  an  orthogonal  matrix  E(z,t). 
Thus , we  have 


(7) 


Ef 


Oz 


row-reduced  with  U being  a k x M matrix  of  rank  k.  (Ever  since  Levin 
(1957),  it  has  been  common  in  the  singular  perturbations  literature  to  use 
block  diagonalization  in  such  situations.  It  seems  more  convenient,  how- 
ever, to  instead  reduce  to  triangular  form.  For  purposes  of  efficient 
computation,  we  note  that  it  is  important  to  use  iterative  procedures  to 
compute  E (cf.  Golub  (1965)).  Defining 


(8) 


A = (0  IM_k)E  and  B = (Ik0)E, 


the  matrix  A(z,t)  annihilates  fQz  and  B(z,t)  is  its  orthogonal 
complement,  i.e., 


(9) 


Af 


Oz 


0,  A3'  = 0,  AA'  = I 


M-k» 


SB' 


Ik  and 


A ' A + B'B  = I 


M* 


Such  matrices  are  convenient  since  they  allow  us  to  represent  any  M vecto 
w in  the  form  w * A'u  + B'v  and  separately  examine  u = Aw  and  v = Ew . 
(We  note  that  Clasen  et  al.  (1976)  used  such  matrices  to  develop  numerical 
schemes  for  solving  stiff  differential  equations.)  It  Is  important  to 


realize  that  our  assumptions  on  f 


imply  that  the  matrix 


Oz 


(10)  S = I3f  B' 

u z 

is  stable,  hence  nonsingular,  because 

Bf0zB'  BB0zA “ 
0 0 


has  k eigenvalues  with  negative  real  parts.  The  stability  condition  (10) 
will  be  especially  critical  in  obtaining  the  boundary  layer  behavior 
needed  at  t - 0. 

We  will  completely  solve  several  large  classes  of  problems,  but  the 
original  problem  must  be  studied  further.  We  shall  show  how  the  outer 
solution  and  the  boundary  layer  correction  can  be  generally  obtained,  but 
our  matching  procedure  for  determining  Z(0,£)  is  not  yet  satisfactory. 

The  Outer  Solution 

Since  the  boundary  layer  correction  n is  asymptotically  negligible 
for  t > 0,  the  representation  (5)  implies  that  the  outer  solution 


(11)  Z(t,€)  -v  l Z , (t  )€ 

J = 0 J 

must  asymptotically  satisfy  the  system 


i 


(12) 


oo 

€ it  = f(Z,t,€)  l f.(Z,t)6‘ 
J = 0 J 
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for  t > 0.  When  £ = 0,  then,  we  have  the  reduced  system  onl — 


(13) 


f0(z0,t)  = 0. 


This  fails  to  determine  ZQ  because  f is  singular.  However,  (9)  and 
(13)  imply  that 


(1H) 


B(Z0,t)f0(A'(Z0,t)A(Z0,t)Z0  + B*(Z0,t)B(Z0,t)Z0,t)  = 


= 0, 


and  since  the  Jacobian  S with  respect  to  BZQ  is  nonsingular,  we  can 
determine  B(Zg,t)ZQ  in  terms  of  the  many  remaining  variables.  We  shall 
assume  that  the  system  (13)  is  consistent.  (An  example  where  the  reduced 


□ □ 


system  is  inconsistent  is  provided  by  fcz^  = -z^  + b,  €z2  = b^,  b2  / 0. 
Then,  the  solution  is  0(i).) 

Vc  shall  obtain  differential  equations  for  ZQ  and  later  terms  Zj 
in  a termwise  fashion.  Rewriting  (12),  we  have 

(15)  = B'(Z0,t)B(Z0,t)f0z(Z0,t)(Z  - ZQ) 

+ (f(Z,t,€)  - fQz(Z0,t)(Z  - ZQ)  ) , 

since  f.  = B ' Bf _ by  (9).  Dif ferentiating  (12)  with  respect  to  t, 

uz  uz 

however,  provides  the  Lagrangian  derivative 


fz(z,t,e)z  + ft(z,t,6)  - ez. 

Thus,  multiplication  by  B implies  that 

(16)  B(Z0,t)f0z(ZQ,t)Z  = -B(Z0,t)ft(Z,t,€)  + B(ZQ,t)(f0z(Z0,t) 

•• 

- fz(z,t,e))z  + tB(z0,t)z. 

Using  (15),  (16),  and  the  invertibility  of  S,  then,  implies  that 

(17)  ^B(Z0,t)f0z(Z0,t)(Z  - ZQ) 

- S"1(Z0,t)B(Z0,t){-ffc(Z,t,6)  + ez  + (f0z(Z0,t)  - fz(Z,t,€))Z 

“ ^Oz(V)(f(Z>t»e)  ' r0z(Z0’t)(Z  ■ Z0))}* 

Finally  substituting  into  (15),  we  have 

(18)  i - A(Z0,t)[-ft(Z,t,€)  + £Z  + (f0z(Z0,t)  - fz(Z,t,6))Z] 

+ g8(Z0,t)[f(Z,t,6)  - f0z(ZQ,t)(Z  - Z0)] 


where 


(19) 


I 


A(Z0,t)  = B,(Z0,t)S"1(Z0,t)B(Z0,t) 
and 

8(Z0,t)  = I - A(Z0,t)f0z(Z0,t). 

? 

We  note  that  the  projection  8=8  satisfies  8B*  = 0 and  Bf„  B = 0. 

Oz 

When  £ = 0,  then,  we  find  that  Zg  must  satisfy  the  nonlinear 
equation 


(20) 


Zg(t)  = - A(Zg,t)f0t(Z0,t)  + B(Z0,t)f1(Zg,t) 


Likewise,  the  coefficient  of  £ in  ( 1 8 ) implies  that  Z^  must  satisfy 
the  nonlinear  equation 


(21) 


Zj(t)  = A(Z0>t)[-fotz(Z0,t)Z1  - flt(Z0,t)  + Z0 


- <WZ0’t)Zl  + flz<Z0>t))Z0] 

+ B(z0,t)t|f0z2(z0,t)zz  + rl!,(z0,t)z1  + f2(z0,t)] 


where  the  quadratic  term  must  be  interpreted  appropriately.  Succeeding 
terms  Zj  will  satisfy  linear  equations  obtained  as  the  coefficient  of 
6^  in  (18).  Using  the.  resulting  differential  equations  (20),  (21),  etc., 
it  follows  that  Zq,  Z^,  ...  will  be  uniquely  obtained  on  some  subinterval 
of  t >_  0 once  their  initial  values  are  specified.  These  initial  values 
by  (12),  must  necessarily  satisfy 

(22)  €fc(0,€)  = f(Z(0,€),0,6). 

Since  the  prescribed  initial  value  z^(0)  will  generally  not  satisfy  (22) 
at  € = 0,  the  need  for  an  initial  boundary  layer  is  obvious. 

Expanding  (22)  and  introducing  1^  = A'A  + B’B,  we  have 

(23)  f0(Z0(0),0)  = 0 

and 

(24)  B0(Z(0,6)  - ZQ(0) ) = - Sg1B0[fgzA0A^(Z(0,e)  - Zg(0) ) 

♦ (f(Z(0,6),0,6)  - f°z(Z(0,6)  - Z0(0))  - €Z(0,6)] 


where 


(D0,S0,f0z,A0)  “ (B,S,f0z?A)|  (ZQ(0)  ,0)  * 

Equatinc  coefficients  of  ^ , J > 0,  in  (24)  implies  that 

B0Zj(°)  - VBofSzVozj(0)  + a j _x 

where  aj_i  is  known  in  terms  of  the  initial  values  Z^O)  of  earlier 
terms.  Thus,  (9)  and  (19)  imply  that 


(26)  Zj(0)  = B0A’A0Zj  (0)  + BqCIj  , J > 0, 
where  BQ  = B(ZQ(0),0),  so  only  the  initial  values 

(27)  VoVj(0) 


need  be  specified  for  each  j > 0.  The  equation  (23)  for  ZQ(0)  also  re- 
mains. (Equation  (14)  implies  that  BQZ0(0)  can  be  uniquely  obtained  in 
terms  of  SQ,  BQ , and  A’,AQZ0(0),  but  the  relationship  will  generally 
be  much  more  complicated  than  the  linear  relationship  (26).) 

The  Boundary  Layer  Correction  and  Matching 
Since  the  outer  solution  Z satisfies  (12),  (5)  implies  that  the 
boundary  layer  correction  II  ( x ,€ ) must  be  a decaying  solution  of  the 
boundary  layer  equation 

(28)  37  = f(Z  + n.CT.e)  - f(Z,6t,€). 


Further,  since  n 0 as  t -*•  °°,  we  must  have 

(29)  n(t,fc)  = -f  [f(Z(€s,€)  + n(s,£),€s,£)  - f (Z(6s  ,6 ) ,€s ,6  ) ]ds  . 

' T 

Before  reconsidering  the  general  problem,  we  consider  some  special  cases. 

a.  The  quasllinear  problem 

We'll  consider  the  quasllinear  problem  with 


(30)  r0(z,t)  - F0(t)z  + G0(t) 

(cf.,  also,  Butuzov  and  Vasil'eva  (1970)).  Here,  fgz,  s» 


A,  and  B are  all  Independent  of  z, 


so 


(31)  Z(t,fc)  = A'(t)A(t)Z(t,€)  + B'(t)B(t)Z(t,fc) 

(cf.  (9)j  noting  that  A and  B are  now  completely  specified).  Further, 
f0(ZQ,t)  = 0 implies  that 

(32)  B(t)Z0(t)  = -S“1(t)[B(t)F0(t)A’(t)A(t)Z0(t)  + GQ (t  ) ] 

and  there  remains  a nonlinear  initial  value  problem  for  A(t)Zg(t)  (cf. 
(31)  and  (20)).  For  higher  order  terms  Zj,  we’ll  now  obtain  BZj  as  a 
linear  function  of  AZj  and  a linear  differential  equation  for  AZj  will 
remain  with  nonhomogeneous  terms  being  functions  of  the  preceding  Z^'s. 
Here,  the  boundary  layer  correction  II  will  satisfy 

(33)  |f  - f0(£t)ii  + { f (z  + n,€T,6)  - f(z,€x,e)  - f0(z  + n,€x) 

+ f0(Z,6x)}. 

Thus,  when  t = 0,  (9)  implies  that 

dno 

(3*0  -^r  = F-0(0)(A’(0)A(0)  + B‘(0)B(0))no 

and  A(0)Fq(0)  = A(0)f0  (0,0)  = 0 Implies  that 

(35)  A(0)nQ(x)  = 0 

since  II q -*  0 as  x *►  ».  From  (3*0,  it  follows  that 

-pj:(B(0)no)  = S(0)B(0)n0, 

hence 

B(0)n0(x)  - eS(0)TB(0)n0(0) 
and,  by  (5),  (9),  and  (35), 

(36)  n0(x)  = B'(0)eS(0)TB(0)(z°(0)  - ZQ(0)). 


w 


We  note  that  (35)  implies  that 


(37)  A(0)Zq(0)  = A(0)z°(0), 

so  we  can  now  integrate  the  initial  value  problem  for  AZQ.  Then,  (32) 
will  provide  BZq  which  by  (31)  and  (36)  imply  the  leading  terms  ZQ 
and  J1  q of  both  the  outer  solution  and  the  boundary  layer  correction. 
We  note,  in  particular,  that  our  formulas  imply  that 

(38)  ZQ(0)  = B(0)A' ( 0) A ( 0)z° ( 0 ) - B ' ( 0 )S-1 ( 0 )GQ ( 0 ) . 


Higher  order  terms  follow  analogously.  Thus,  we  have 


dn, 

-ar  " V0)ni + 


where  y^(t)  is  known  and  exponentially  decaying.  Clearly,  A(0)Fq(0)  = 0 
implies  that 


(39) 


A(0)n1(x)  - A ( 0 ) 


(s  )ds 


and  this  provides  the  initial  value  A(0)Z^(0)  = A(0)(zq1  - 11^(0))  needed 
to  completely  determine  A(t)Z^(t)  and  Z^(t).  Continuing,  we  have 


(^0)  B ( 0 ) II 1 ( x ) = eS(0)TB(0)(z01  - Zx(0)) 

. T *oo 

+ J eS(0) (t"s)B(0)[F0(0)A' (O)A(O) J Y1(r)dr  + Y1(s)]ds 

and 


(Jji)  n1(i)  = A,(o)A(o)n1(x)  + B,(o)B(o)n1(x) 

is  completely  specified.  In  this  way , all  terms  of  the  asymptotic  solution 
to  the  quasilinear  problem  (30)  are  completely  and  uniquely  determined. 

b.  The  classic  problem 


A much-studied  nonlinear  initial  value  problem 


(42)' 


^clt  = v(x>y»t>e)»  y(o,e)  = y°(£) 

^ = u(x,y,t,6),  x(0,6)  = x°(fc) 


can  be  analyzed  under  the  principal  assumption  that  the  reduced  problem 


(43) 


10  = vQ (x ,y  j t ) 

|^r  = uQ(x ,y  ,t ) , x (0)  = x° ( 0 ) 


has  a unique  solution  (XQ(t),  YQ(t))  along  which  v0^(XQ(t),  YQ(t),t) 
< 0 for  t > 0 (cf . , e.g.,  O'Malley  (1974)).  For  simplicity,  we'll  as- 
sume that  x and  y are  scalars  and  that 


(44) 


vQy(x,y,t)  < 0 


everywhere.  The  problem  (42)  can  be  put  into  the  form  (l)-(2)  by  intro- 
ducing z = , so  that  we  have 


(45) 


£z  = f (z,t,€) 


v(x,y,t,6) 
€u(x,y ,t,6) 


z (0,6) 


0 

z 


(6). 


The  limiting  system  f(ZQ,t,0)  = 0 reduces  to 
(46)  v0(XQ,Y0,t)  = 0 


and  the  Jacobian  matrix 


(47) 


has  rank  one  under  hypothesis  (44).  (Here 
of  eigenvectors,  since  the  null  space  of 
one-dimensional. ) 

Because  fQz  is  already  row-reduced , 
A = (0  1),  and  B = (1  0),  so  that  y 

problem  (46)  can  be  solved  for  YQ  = <KXQ, 


, there  cannot  be  a deficiency 
f q and  its  transpose  are  both 

we  can  take  E = { q ^}  , 

■ Bz  and  x = Az.  The  reduced 
t)  since  v0y  is  nonsingular. 


but  (46)  fails  to  determine  XQ.  The  critical  assumption  (Ij4)  corresponds 
to  the  earlier  hypothesis  (10),  since  S = vQy  < 0. 

For  this  problem,  the  leading  term  Zg  of  the  outer  solution  satis- 
fies the  equations 


(48) 


v0xu0>- 


U 


0 


(a  specialization  of  (20)),  Indeed,  we  can  rewrite  this  system  as 

d * 

^-v0(X0,Y0,t)  = 0 with  XQ  = Uq,  so  the  initial  condition  vq(Xq(0), 
Yq(0),0)  = 0 reduces  the  system  to  the  familiar  reduced  problem  (43) 
provided  Xg(0)  = x^(0). 

If  we  further  set  the  boundary  layer  correction  II  = |^|,  the  limitin 
form  of  the  boundary  layer  system  (28)  reduces  to  the  pair  of  equations 


(49) 


dn0 

dF  - W0' 


W0)  + n0’0^  * dT  °' 


Moreover,  the  decay  required  as  x -*•  °°  implies  that  £n(x)  = 0,  so  we 

0 u 

must  have  Xg(0)  = x (0)  and  the  limiting  outer  solution  (XQ(t ) ,YQ(t ) ) 
is  uniquely  specified.  Thus,  there  remains  the  problem 

(5°)  = vo(x°(O),0(x°(O),O)  + n0,0),  nQ(0)  = y°(0)  - $(x°(0),0) 

The  Dini-Hukahara  theorem  (cf.  Hille  (1969))  and  the  negativity  of  v0y» 
imply  the  existence  of  a unique  decaying  solution  ng(t).  Thus,  the 
leading  boundary  layer  correction  term  nQ  is  completely  determined. 
Succeeding  terms  follow  i a straight  forward  manner. 

We  note  that  a considerably  more  complicated  vector  problem  would 
result  if  vQy  had  zero  eigenvalues. 

c.  The  original  problem 

For  the  problem  (1)— (2) , the  Initial  value  BqZq(0)  of  the  BQ 
portion  of  the  leading  term  of  the  outer  solution  Is  restricted  by  (14). 

It  restricts  the  leading  boundary  layer  term  nQ  since  it  must  satisfy 
Bgn0(0)  = Eq(z°(0)  - Zq(0))  (cf.  (5)).  Alternatively  the  initial  value 
AgZg(O)  is  free  as  far  as  the  outer  solution  is  concerned,  though  it  would 
be  determined  through  AqIIq(0)  = Ag(z^(0)  - ZQ(0)). 

Multiplying  (28)  by  BQ  and  integrating  from  0 to  x implies  that 

SnT  « f t Sg(x**s) 

B^B0ri(x,O  = B£e  U (zU(S)  - Z(0,€))  + BJJ  e BQ 


(51) 


[f(Z(ts,6)  + A^Af)n(s,€)  + BJ)Bon(3,e),eS,C) 
- f(Z(es1€),€s,6)  - f°zB'B0n(3,£)]ds 


while  (29)  Implies  that 


(52) 


Pi'QkQ\[(T  ,£)  ~HQAQ 


CT(Z(6r,t)  + 


AoAon(r,e) 


B’Bon(r,£) 


£r,fc)  - f('Z(€r,6),6r,6)]dr. 


Putting  these  together,  we  get  the  new  integral  equation 


(53)  n(x,€ ) = ji°(t)  + F(n) 

where 


n (t)  = 


S0T 

B;e  Bq(2 


(t)  - Z(0,6)) 


and 


F(n)  = B 


Sq(t-s) 
e 


B0[f(Z(6s,e) 

- r^BjB0n(s,6)]d3  - aja f 

J T 


+ II  (s,<=  ) ,fcs  ,€  ) - f (z(6s  ) ,€s  ,€  ) 

[f(Z(€r,6)  + n(r,e),€r,€) 


- f(Z(6r,6),fcr,fc)]dr. 


We  recall  that  for  any  initial  value  Z(0,£),  consistent  with  (12), 
Z(t,£)  is  completely  determined.  Likewise,  the  exponential  decay  cf 

v 

e implies  that  for  any  arbitrary  initial  vector  BqTT  ( 0 , 6)  a k-di- 
mensional  manifold  of  exponentially  decaying  solutions  of  (53)  is  ob- 
tainable by  successive  approximations.  We  simply  define 


(5*0 


rjl  (J+1)(T  ,t)  = n °(t  ) + F(JI  (J)), 

Sand 

n (t  ,e)  = lim  n (j)(t  ,t) 

>°° 


3 > 0 


(cf.  Erdelyi  (I96U) ) . As  for  the  quasilinear  and  classic  problems,  we 
are  prescribing  the  vector  BqBJ1(0,£)  (of  rank  <_  k)  and  letting  the 


vector  A-AgnCO.C)  (of  rank  <_  M - k)  be  determined  by  the  exponentially 
decaying  solution  constructed  (cf.  (51)  and  (52)).  Vasil'eva  (1975) 
expresses  this  relationship  by  writing  an  expression  like  AqAqII(0,€) 

= «(B^Bon  (0,c)  ,Z(0,<i) ) . This  actually  relates  initial  values  of  the 
outer  solution  s3nce  it  implies  that 

(55)  A£A0(s°(t)  - Z(0f6))  = *(B£B0(z°(€)  - Z(0,€)),  Z(0,6)). 


(56) 


It  is  instructive  to  further  consider  the  asymptotic  expansion 

on 


n( t ,£)  'v  l n,  (t)€i 

j=o  J 


for  the  boundary  layer  correction.  Clearly,  n q must  satisfy  the 
nonlinear  system 


dH 


(57) 


a r = W°>  + "o>0) 


and  we  have 
(58) 


0 rT  Sq(t-s) 

n0(T)  = nu(x ) + j B^e 

" f0zB0B0n0(s) "ds 


B0[f0(z0(0)  + n0(s),0) 


.00 


AoAoj 


f0(Z0(0)  + no(r),0)dr. 


We  note,  in  particular,  that  we'll  generally  have  Ao^O^T^  ^ unlike 
for  the  quasilinear  or  classic  problems  where  aq**Oz  = 0*  Moreover,  II ^ 
will  satisfy  a linear  system 


(59) 


dn, 

= f0z(Z0(0)  + no,0)ITl  + a0(x)II0 


where 


°o(T)no 


lim  {-  |{f(z  + n0,£x,e)  - f(zQ(o)  + n0,o,o) 


f(Z,€x,6)}. 


Since  the  coefficient  matrix  CQz  has  k stable  eigenvalues,  we 
can  anticipate  a k-dinensional  manifold  of  initial  values  1^(0)  cor- 
responding to  decaying  solutions  as  x This  allows  us  to  prescribe 

BqB0 n i ( 0 ) while  the  remaining  components  AqAqN^O)  of  the  initial 


vector  become  specified,  (We  recall  that  Vasil 'eva  (1975)  seems  to  re- 
late these  initial  values  in  terms  or  the  M - k eigenvectors  cpanninr 

the  null  space  of  f ' . ) Analogous  considerations  apply  for  higher  order 

u z 

terms . 

A further  problem  remains;  namely,  the  initial  value  Z(0,£)  of 
the  outer  solution  is  still  imbedded  in  the  nonlinear  equations  for  the 
outer  solution  and  for  the  boundary  layer  correction  (cf.  (12)  and  (14) 
with  t = 0 and  (55)).  We  shall,  for  the  time  being,  simply  assume  that 
the  resulting  nonlinear  equations  have  a unique  solution  Z(0,fc)  as  they 
do  in  the  special  cases  already  considered  (cf.  also  Vasil'eva  (1975) ) - 
More  critically,  it  is  important  to  develop  numerical  methods  for  stiff 
systems  based  cn  this  analysis.  Such  work  is  in  progress. 
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